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Abstract 

Based on the UrQMD transport model, the transverse momentum and the rapidity dependence 
of the Hanbury-Brown-Twiss (HBT) radii Rl, Ro, Rs as well as the cross term Rol at SPS ener- 
gies are investigated and compared with the experimental NA49 and CERES data. The rapidity 
dependence of the Rl, Ro, Rs is weak while the Rol is significantly increased at large rapidities 
and small transverse momenta. The HBT "life-time" issue (the phenomenon that the calculated 



Rq — R? s value is larger than the correspondingly extracted experimental data) is also present 
at SPS energies. 
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The high temperature phase of Quantum Chromodynamics (QCD), i.e., the Quark Gluon 
Plasma (QGP) has been investigated theoretically and experimentally for many years with 
energies in y/s ranging from around 2 GeV up to 200 GeV. Since the (phase) transition from 
hadrons to quarks might take place only in a small sub-volume and for a short timespan 
compared to the overall process of the heavy ion collisions (HICs), it is an intricate but 
essential task to explore the space-time structure of the region of homogeneity in HICs. For- 
tunately, a well-established technique called Hanbury-Brown-Twiss interferometry (HBT) 
has been developed, which is now extensively used in the heavy-ion community to 
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the particle emission source for HICs with ener gies 
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35]. Studies on various species of two-particles other 



than identical charged pions have been published or are being pursued (see, e.g. Refs. 
and the references therein) . After more than two decades of progress aiming at relativistic 
HICs, many basic and important systematic features have been discovered by femtoscopic 
measurements at SIS, AGS, SPS, and RHIC 0, 16, 23, 24]: the dependence on system size, 



collision centrality, rapidity, transverse momentum and particle mass (or species). 



In our recent work 



361 ]. we have investigated theoretically the system size- and the cen- 



trality dependence of the HBT radius parameters Rl, Ro, and Rs at RHIC energies. It was 
found that the calculations are well in line with experimental data, except that the predicted 
Ro values in central collisions (besides the elementary p+p collisions) are slightly larger than 
the experimental data. Thus, the R Q j R s ratio is larger than that extracted from data. This 
phenomenon is sometimes dubbed the "HBT-puzzle" and reported by most of the other 
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37J, lag). In Ref. |4l| A Multi-Phase Transport model (AMPT) provides 



a good fit to data after including a partonic cascade with a certain parton-scattering cross 
section, while in Ref. 14j the Hadronic Rescattering Model (HRM), which is only based 



on hadronic rescattering with instantaneous collisions, can interpret HBT observables as 
well. Furthermore, J. G. Cramer et. al. 39] presented one relativisitc quantum-mechanical 
treatment of refractive and opacity effects on two-pion correlations and found a chiral phase 
transition is required in the dense medium in order to solve the HBT "puzzle" at RHIC, 
which was further supported b y S . Pratt in 40J by considering interactions with mean field 
in classical trajectories. Ref. (33] argued that the origin of this "HBT- puzzle" might be 
multifaceted. In this work, we focus our investigation on HICs at lower energies — the 
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SPS energy regime — partly in order to better understand the origin of the large R Q j R s in 
present models. Furthermore, in this energy region, recently the rapidity dependence of the 
HBT-radii of the pion source was investigated systematically by two experimental collabo- 



rations NA49 



35, 



421 and CERES 



It is expected that the rapidity (or pseudorapidity) 
dependence of the HBT parameters might provide essential information on the localization 
of hot and dense nuclear matter and hence it is worthwhile to investigate theoretically. 

By using the Ultra-relativistic Quantum Molecular Dynamics (UrQMD, v2.2) transport 
model (employin g ha dronic and string degrees of freedom) (for details, the reader is referred 



to Refs. L43L 



44 
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461 ]) and the analyzing program CRAB (v3.0/5) 



47|, |48|, |49[ as tools, 



we compare the calculated transverse momentum and rapidity dependence of HBT-radii 
with the (preliminary) experimental data at SPS energies Q, 42]. The HBT-radii are 
represented by using the Bertsch-Pratt three-dimensional convention, which is also called 
the longitudinal co-moving system or "Out-Side-Long" system. In this work, the three- 
dimensional correlation function is fit with the standard Gaussian form: 



C(qo, qs, Ql) = 1 + \exp(-R 2 L q 2 L - R 2 q 2 - R 2 s q 2 s - 2R 2 OL q q L ), (1) 

in which g« and Ri are the components of the pair momentum difference q = P2 — Pi and the 
homogeneity length (the HBT-radii) in the i direction, respectively. The A parameter is the 
incoherence factor which lies for Bose-Einstein statistics between (complete coherence) and 
1 (complete incoherence). Experimentally, the A parameter is affected by many factors e.g. 
non-pionic contaminations, (long-lived) resonances, and details of the Coulomb correction. 
Theoretically, the A parameter is mainly suppressed by the presence of mediate/long-lived 
resonances since pions from the decay of these resonances come from a larger source than the 
directly emitted pions, which results in an overall suppression of the correlation function. 
The Rql represents the cross-term, which is supposed to vanish at mid-rapidity for symmet- 



ric systems. At large rapidities, Rol deviates from zero [50j, |5l|], which will be investigated 
in this paper as well. In the present UrQMD calculations the potential interactions are not 
considered (the cascade mode) due to the oppressive computing time. Coulomb final state 
interactions (FSI) are not taken into account in the program CRAB and the fitting process. 

In Figs. [Hand [2] we show the transverse momentum kx dependence (ky = (p 1T + p 2T )/2) 
of the HBT- radii Rl (top plots), R<j (middle plots), and Rs (bottom plots) in Pb+Pb 
reactions at beam energies E b = 20, 30, 40, 80, and 160A GeV, and in Pb+Au reactions at 
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FIG. 1: (Color Online) Transverse momentum dependence of the Pratt-radii in central Pb+Pb 
collisions at = 20, 30, 40, 80, and 160A GeV and mid-rapidity < Y n7T < 0.5. Preliminary 



NA49 data are taken from [42j. At Ef, = 160A GeV the calculation results of the HBT-radii with 
a vanishing formation time for strings are also presented. As a result, the /cy-dependence of the 
HBT radii becomes steeper, and the values of Rs are increased and approach the calculated Ro 
values. 

beam energies Ej, = 40, 80, and 160A GeV (plots from left to right) , respectively. In Fig. [1] 
the results are compared with preliminary NA49 data [42] (< 7.2% of the total cross section 
a T ). The pair-rapidity < Y T7T < 0.5 (Y W7T = |log( ^^^|j|^jp is the two-pion rapidity 
with energies E\ and E 2 and longitudinal momenta p\\% and pp in the center-of-mass system) 
is chosen for all reactions. In Fig. [2], the results are compared with CERES data 18j (the 
< 5% most central collisions). The \Y nn \ < 0.25 for 40A GeV, -0.5 < Y n7T < for 80A GeV, 
and —1.0 < Y nn < —0.5 for 160A GeV, respectively. Furthermore, the n~ — n~ correlations 
are calculated in Fig. [T] while the two-charged-pion correlations (including two-7r~ and two- 
7r + mesons) are calculated in Fig. [2] (same as experimental outputs). The statistical errors 
are shown as well. The calculated HBT radii at 160 A GeV by adopting zero formation time 
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FIG. 2: (Color Online) Transverse momentum dependence of the Pratt-radii in central Pb+Au 
collisions at E b = 40 flY^I < 0.25), 80 (-0.5 < Y n7T < 0), and 160A GeV (-1.0 < Y n7T < -0.5). 



CERES data are taken from 



lfil- 



for strings (r s = fm/c) will be explained together with Fig. [TTJ Firstly, it is very interesting 
to see that the present calculations can reproduce the /c^-dependence of HBT radii Rl and 
Rs fairly well shown in Figs. [I] and [2j Only at small &t values, the calculated Rl and 
Rs values are seen up to 25% lower than data. While for the Ro values, the calculations 
are shown larger than both NA49 and CERES data especially at relatively large kx- This 
deviation appears strongest in most central collisions, which was also seen at RHIC energies 



3fj. By comparing the NA49 data with the CERES data for central Pb+Pb (a/a T < 7.2%) 
and central Pb+Au (<j/<jt < 5%) collisions, one observes that the CERES Ro data are 
somewhat smaller than the recent published NA49 data [42j especially at large fcy and low 
beam energy Ej, = 40A GeV although the recent NA49 data at large fcr have already been 
driven down visibly when comparing with those preliminary data in the previous publication 



35[. The difference between NA49 and CERES data is obviously not due to the differently 
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FIG. 3: (Color Online) /cr-dependence of the cross term Rq L . Central Pb+Au collisions at Ef, = 40 
(\Y nn \ < 0.25), 80 (-0.5 < Y nn < 0), and 160A GeV (-1.0 < Y wn < -0.5) are chosen for comparing 
with the CERES data taken from [lj. 



charged pion species since, similar to the HBT results at RHIC |28l . |30j , we also find that 
the difference between the HBT-radii of two positively charged and two negatively charged 
pions is quite small at SPS energies. The origin of these differences between data are still 
not quite clear. 

In addition to the difference between both experimental data in centrality, the selected 
rapidity cuts are also different in both experiments. The effect of the rapidity-cut on the 
hr- dependence of the cross term Rq L is shown in Fig. [3] in comparison with the CERES 
data. In line with the data, an increase of the absolute value of Rql ^ s seen a ^ small kx 
with the increase of beam energy. This increase is mainly due to the shift of the rapidity 
cut away from mid-rapidity with the increase of the beam energy. This effect will also be 
shown in Figs. [7] and [HI When kr approaches MeV/c, the value of the cross term Rq L also 



approaches zero, as one expects [511 ]. 



Preliminary results on the rapidity dependence of n n Bose-Einstein correlations in 
central Pb+Pb collisions at 20A, 30A, A0A, 80A, and 160A GeV measured by the NA49 
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collaboration are available 42j • In Figs. H] and [5] we show the rapidity dependence of the 
HBT radii Rl (upper plots), Ro (middle plots), and Rs (bottom plots) within two transverse 
momentum bins: < < lOOMeV/c and lOOMeV/c < kr < 200MeV/c, respectively. Here 
as well as below we only choose three cases with beam energies 20, 40, and 160A GeV as 
examples. For Rs values shown in Figs. H] and the Rs of NA49 data decrease with 
increasing rapidity slowly which is also observed in our calculations despite slightly smaller 
radii. For Ro results it is seen that the rapidity dependence is weaker than for Rs- The 
comparison of our calculated results with NA49 data for the rapidity dependence of the Ro 
values is reasonably well, with the deviation up to 20%. For Rl results, we find a stronger 
rapidity dependence in the low transverse momentum bin < kr < lOOMeV/c while it is 
reduced at large momenta. At transverse momenta < kr < lOOMeV/c and at low beam 
energy Ej, = 20 A GeV, the Rl value first drops down and then rises again with the increase 
of rapidity, which is observed as well in our calculations. While at large beam energy 160A 
GeV a different calculated rapidity dependence of Rl values compared with the preliminary 
experimental data is found, that is, with the increase of rapidity the calculated Rl value 
first rises and then drops while the trend of experimental data is still preserved similar to 
the reactions at lower beam energies although the minimum shifts to a larger rapidity. 

As we know, in UrQMD transport model particles are produced via the decay of a meson- 
or baryon-resonance or by color-strings excitation and fragmentation (the antibaryon-baryon 
annihilation yields negligible contributions). Keep in mind that the two ways of particle 
production are different: the former is based on the quasi-particle level (with respect to 
the resonance masses m < 2.5 GeV), while the latter happens at large two-particle center- 
of-mass energies ( yj $nn ^ 2 GeV). Therefore, the string process dominates at the early 
stage of HICs with large beam energies. While the decay of resonance takes place over the 
whole evolution of the system. In order to understand the dissimilar rapidity-dependence of 
Rl at large beam energies, we show in Fig. [6] the single-pion rapidity (Y cm = |log( ^^||| )) 
distribution of the ratio of the pion numbers at freeze-out from the decay of resonances and 
strings. Firstly, we see that the mean ratios have exceeded unit which means the decay of 
resonances dominates at freeze-out and is also of importance to the HBT parameters. In 
this figure we also find an interesting beam-energy dependence of the ratio as a function of 
rapidity, e.g., at Eb = 20A GeV, it first drops from mid-rapidity to projectile-target rapidity 
region and then rises further with increasing rapidity. While at the largest energy point of 
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FIG. 4: (Color Online) Rapidity dependence of the HBT-radii Rl (upper plots), Ro (middle plots), 
and Rs (bottom plots) within transverse momentum bin < hr < lOOMeV/c at beam energies 20, 
40, and 160A GeV. The calculations are shown with lines. Preliminary experimental NA49 data 
are taken from 42( (stars). 



SPS, Ef, = 160A GeV, the ratio at mid-rapidity is suppressed due that much more pions 
produced by string fragmentation are treated to transversely emit with small longitudinal 
velocity. It is known that the case when more pions emitted via the decay of resonance leads 
to the postpone of pion freeze-out and hence extend the volume of source. While larger 
fraction of pions via the decay of strings means the freeze-out takes place at earlier stage 
of overall process and hence a smaller volume of source. Therefore, the cancelation of the 
two mechanisms of pion freeze-out determines the explicit rapidity dependence of the HBT 
radii. 

Next, we show the rapidity dependence of Rol at transverse momenta < < 
lOOMeV/c (in Fig.[?D and lOOMeV/c < k T < 200MeV/c (in Fig. E>, in Pb+Pb reactions at 20 
(top), 40 (middle), and 160A GeV (bottom). It is seen clearly that the cross term increases 
with the increase of rapidity which is expected since the longitudinally boost-invariance is 
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FIG. 5: (Color Online) Rapidity dependence of the HBT-radii Rl (upper plots), Ro (middle plots), 
and Rs (bottom plots) within transverse momentum bin 100 < kx < 200MeV/c at beam energies 
20, 40, and 160A GeV. The calculations are shown with lines. Preliminary experimental NA49 
data are taken from 42| (stars). 



broken in the forward and backward rapidity regions [5JJ. We find in our calculations that 
the inclusion of this cross term modifies the HBT radii Rl, Ro, and Rs only weakly when 
adopting the rapidity cut < Y nn < 0.5. However, it is necessary to include this cross term 
at large rapidities and small &t because then it deviates from zero, as shown in Figs. El d 
and [HJ One also observes that the present calculations reproduce the rapidity dependence 
of Rol values reasonably well in the rapidity region < < 2. 

Figs. [9] and [10] show the rapidity dependence of the quantity a/ R q — R s in central Pb+Pb 
reactions at < kx < lOOMeV/c and 100 < kx < 200MeV/c, respectively. In certain 



scenarios 52], the so-called "duration time" At of the pion source is proportional to this 
quantity at a fixed transverse momentum and has been extensively discussed in the HBT 
community. However, one should be careful, when interpreting this quantity as "duration 



time" because this identification only holds in the absence of flow and opacity effects 52]. 
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FIG. 6: (Color Online) Single-particle rapidity (Y cm ) distribution of the ratio of the pions at freeze- 
out from the decay of resonances and strings. Three cases at beam energies 20, 40, and 160A GeV 
are chosen. A similar beam-energy dependence as Fig. [4] in Rl is seen, (see text for details) 



Mainly due to a smaller calculated Rs in both transverse momentum bins, the calculated 



a/ Rq — R 2 S values are larger than the extracted experimental data at the rapidity interval 
considered. It is interesting to see that the calculation seems to approach the data at large 
rapidities Y nn = 1.5 ~ 2.0. In order to solve the HBT "puzzle", the treatment of the 
interactions during the whole process of HICs should be important. One should notice that 
the largest deviation of the calculated HBT radii from experimental data comes from the 
pion correlation at mid-rapidities in the most central collisions, where the interactions are 
strongest. Owing to the omission of treatment on string-string interactions in UrQMD at 
early stage, one can argue that the pressure in the early stage might be too small, resulting in 
a delayed source break-up. The small elliptic flow calculated at RHIC energies is attributed 
to the same origin 46(. The mean free path of partons (or formation time of particles) from 
the hot midrapidity region was observed closely linked to the strength of the elliptic flow 53] 
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FIG. 7: (Color Online) Rapidity dependence of Rol calculations (lines) at transverse momenta 
k T < lOOMeV/c in Pb+Pb reactions at 20 (top), 40 (middle), and 160A GeV (bottom), compared 



to NA49 preliminary data 



42j (stars). 



and should also modify the HBT results. In Fig. [T] (in central collisions at Ej, = 160A GeV) 
we show the HBT radii with a zero formation time of particles as well. Due to the reduced 
mean-free-path of partons and hence much larger transverse pressure at early stage, the kx 
fall-off of the HBT-radii is seen even stronger than that found in the data which implies 
stronger coordinate-momentum correlation in the modified calculations. The R$ rises more 
strongly than the other two HBT radii and approaches the Ro value. In Fig.[TT]the extracted 
a/ Rq — R 2 S values from NA49 as well as from CERES data for central collisions as a function 
of at 160A GeV are compared with calculations of central Pb + Pb collisions with the 
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FIG. 8: (Color Online) Rapidity dependence of Rol calculations (lines) at transverse momenta 
lOOMeV/c < k T < 200MeV/c in Pb+Pb reactions at 20 (top), 40 (middle), and 160A GeV 



(bottom), compared to NA49 preliminary data 



4J (stars). 



default and the reduced formation times of strings. Undoubtedly, the reduced formation time 
of particles leads to smaller a/ Rq — R s values. At kx < lOOMeV/c, the calculated R$ values 
are even larger than Rq values. It is noticed that the CERES Rs data at kj- < 500MeV/c 
are also larger than Rs data. It should be noted that, although a shorter formation time 
is apt to explain the so-called "HBT-puzzle" , as well the elliptic flow, the absolute values 
of HBT-radii are not in line with the data (shown in Fig. [[]). A more careful treatment on 
string-string and/or hadron- string interactions as well as the decay of resonances are essential 
to consistently explain both the single-particle flows and the two-particle correlations. 
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FIG. 9: (Color Online) Rapidity dependence of the calculated quantity y R Q — R s (lines) at kx < 
lOOMeV/c with beam energies 20 (top), 40 (middle), and 160A GeV (bottom). The preliminary 



experimental data are illustrated with stars 



42|. 



To summarize, by using the CRAB program, we analyzed the energy and the rapidity 
dependence of the HBT-parameters Rl, Ro, Rs as well as the cross term Rol at SPS energies 
in collisions calculated within the UrQMD transport approach. As a whole, the obtained 
transverse-momentum and rapidity dependences of the HBT-parameters are in reasonable 
agreement with the experimental NA49 and CERES data, except that the calculated Ro at 
large kr and Rs at small kr are somewhat larger than data. The rapidity dependence of 
the HBT radii Rl, Ro, Rs is relatively weak while that of the cross term Rol is significant. 
Comparing our calculations and the experimental data, the so-called "HBT-puzzle" reoccurs 
at SPS energies — the calculated a/ R — R s values at small kx are somewhat larger than the 
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FIG. 10: (Color Online) Rapidity dependence of the calculated quantity y Rq — Rg (lines) at 
100 < k T < 200MeV/c with beam energies 20 (top), 40 (middle), and 160A GeV (bottom). The 
preliminary experimental data are illustrated with stars [42]. 

experimental data. This observation might hint to several open questions in the treatment 
of the dynamics of HICs, such as the omission of the string-string interactions, the omission 
of potential interactions of hadrons during the evolution and the involved treatment of 
resonance production and lifetime in transport models. 
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